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Abstract. The deexcitation of exotic hydrogen atoms in highly excited states in collisions with hydrogen 
molecules has been studied using the classical-trajectory Monte Carlo method. The Coulomb transitions 
with large change of principal quantum number n have been found to be the dominant collisional deex- 
citation mechanism at high n. The molecular structure of the hydrogen target is shown to be essential 
for the dominance of transitions with large An. The external Auger effect has been studied in the eikonal 
approximation. The resulting partial wave cross-sections are consistent with unitarity and provide a more 
reliable input for cascade calculations than the previously used Born approximation. 
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PACS. 34.50.-s Scattering of atoms and molecules 
mesons, muons, and other unusual particles) 



1 Introduction 

Exotic hydrogen atoms x~p {x~ = fi~ ,Tr^ , K~ ,p) are 
formed in highly excited states with the principal quan- 
tum number n \/Jj^yp/rn^ where ^xp is the reduced 
mass of the exotic atom For a long time the ini- 

tial stage of the atomic cascade remained poorly under- 
stood despite a substantial progress in theoretical and ex- 
perimental studies (see |^-|^ and references therein). In 
particular, the dominant collisional deexcitation mecha- 
nism was unclear for 40 years since the so-called chemical 
deexcitation was introduced in j^] as a phenomenological 
solution to the problem of the cascade time at high n (the 
external Auger effect alone would give much longer cas- 
cade times). A shortage of experimental data related to 
the initial stage of the atomic cascade hindered theoret- 
ical studies of this problem. The experimental situation, 
however, changed recently as more data on the atomic 
cascades in exotic hydrogen atoms at low density became 
available. The cascade time of antiprotonic hydrogen mea- 
sured by the OBELIX collaboration in the density 
range 3 — 150 mbar was found to be significantly shorter 
than the prediction of the conventional cascade model . 
The new experimental results on the atomic cascade in 
muonic hydrogen from the PSI experiment Q provided 
detailed information not only on the cascade time, but 
also on the energy distribution at the end of the cascade, 
which at low density is actually preserved from the initial 
stage after the fast radiative deexcitation takes over the 
collisional processes. 

The goal of this paper is to investigate the collisional 
deexcitation mechanisms for highly excited exotic atoms. 



36.10.-k Exotic atoms and molecules (containing 



In particular, we are interested in the role of the Coulomb 
acceleration in highly excited states and in the competi- 
tion between the acceleration and slowing down in quasi- 
elastic collisions. Both molecular and atomic hydrogen tar- 
gets were used in our calculations in order to investigate 
the role of molecular effects. 

The paper is organized as follows. The classical-trajectory 
Monte Carlo method is described in Section |[ The results 
of calculations of Coulomb, Stark, and transport cross- 
sections for the and pp atoms are presented in Sec- 
tion ||. The Auger deexcitation is discussed in Section ||. 
The conclusions are summarized in Section |^. 

Unless otherwise stated, atomic units {h — e — me — 
1) are used throughout this paper. The unit of cross- 



section is Gq ~ 2.8 • 10 cm^, where qq — 
electron Bohr radius. 



is the 



2 Classical-trajectory Monte Carlo calculation 
2.1 Effective potential 

In the beginning of the atomic cascade, where many nlm- 
states are involved in the collisions, classical mechanics 
is expected to be a good approximation. To study the 
scattering of exotic hydrogen atoms from hydrogen atoms 
or molecules we use a classical-trajectory Monte Carlo 
model. The following degrees of freedom are included in 
the model: the constituents of the exotic atom (x~ = 
, TT^ , , p and the proton) and the hydrogen atoms are 
treated as classical particles. The electrons are assumed 



2 



T.S. Jensen, V.E. Markushin: CoUisional deexcitation of exotic hydrogen atoms. I. Cross-sections 




H 



Fig. 1. Coordinates for the three-body subsystem. C is the 
center of mass of the x~p. 
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Fig. 2. Coordinates for the four-body system. 



to have fixed charge distributions corresponding to the Is 
atomic state around the protons in the hydrogen atoms. 

We describe the exotic atom as a classical two-body 
system with the potential 



1 



(1) 



The exotic atom interacts with two hydrogen atoms whose 
electron distributions are assumed to be frozen in the 
ground atomic state (see Figure |^ for notation) : 
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(2) 



The interaction between the hydrogen atoms is described 
by the Morse potential 



I^hh(Rhh) - i?,(e-«(«™-«") - 1) 



(3) 



where = 4.75 eV, a ^ 1.03, and Rq = 1.4 g. The 
effective potential for the p+H+H system (see Figure ^ 
has the form 

V = V,-pir) + V;-p_^(r, Ri) + ^-^^^(r, R2) 
+ Vhh(Rhh). (4) 



2.2 Method of calculation 

The classical equations of motion corresponding to the 
effective potential (W) were solved using a fourth-order 



Runge-Kutta method. The initial conditions were defined 
as follows. Given the initial principal quantum number rii 
and the orbital angular momentum li of the x^p atom, 
the initial classical state was generated as a classical Ke- 
pler orbit with the total CMS energy E^p and the classical 
angular momentum Ic'. 
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(5) 
(6) 



The orbit was oriented randomly in space, and the or- 
bital x~p motion was set at a random time within the 
period. The hydrogen atoms in the target molecule were 
set at the equilibrium distance i?o, and the molecule was 
randomly oriented in space. The impact parameter p of 
the x^p atom was selected with a uniform distribution in 
the interval (0,pniax), as discussed below. The accuracy 
of the numerical calculations was controlled by checking 
the conservation of total energy and angular momentum. 
Instead of requiring convergence for every individual tra- 
jectory, we used the global criteria that the cross-sections 
for the various processes (see below) were stable within the 
statistical errors against further increase in the numerical 
accuracy for each collision. 

The final atomic state was determined when the dis- 
tance between x~p and the hydrogen atoms after the colli- 
sion was larger than lOao. The final atomic state with the 
energy Exp and the angular momentum Ic was identified 
as correspondingto the final njlf state according to the 
rules similar to (|5[|g): 
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(7) 
(8) 



In addition to the quantum numbers n f and I f , the CMS 
scattering angle and the excitation energy of the target 
^-Etargct were also obtained. For the purpose of cascade 
calculations, we are mainly interested in the reaction chan- 
nels that include the x~p atomic states: 

{{x^p)nflj + H2 
{x-p)nflf +^11 (9) 
ix~p)nflf + H + H 

Other possible channels are the breakup reactions 

- +P + K 



(.X P)ni 



H2 



p + H + H (10) 



and the formation of [x H)„^;^ ions 

(x-p)„,;, +H2 ^ {x-ll)nfl, + P + H. (11) 

An example of a collision that results in Coulomb de- 
excitation of the fi~p and dissociation of the hydrogen 
molecule is shown in Figure ^. 

For a given initial n of the x^p and laboratory kinetic 
energy, Ti, a set of the impact parameters pi, i — 1, K 
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Fig. 3. An example of a ii~p + H2 collision with impact 
parameter p = ao resulting in Coulomb deexcitation of the fi~p 
and dissociation of the H2. The exotic atom with laboratory 
kinetic energy Ti = 1 eV, Ui — 14, and h = 13 enters from the 
left, the hydrogen molecule from the right. The trajectories are 
plotted in the CMS. In the final state, the /i~p has n/ = 10, 
If ^7, and T/ = 4.3 eV. 



with a uniform distribution in the interval (0,pmax) was 
generated. The value Pmax = 5 + '^nf/fj,xp was found to be 
suitable for all cases concerned. The cross-sections were 
obtained from the computed set of trajectories using the 
following procedure. Let be the probability that the re- 
action channel a corresponds to the final state in collision 
i: 



pa _ 



1 , if a occurred 
, otherwise 



(12) 



The cross-section for the reaction channel a is given by 



1 



K 



1=1 



(13) 



The differential cross-sections are determined in a similar 
way by binning the corresponding intervals of variables 
like z = cos6', where 9 is the CMS scattering angle, and 
the target excitation energy Z\i?target- For instance, the 
differential cross-section da{z)/dz is calculated using the 
relation 



dcr(z) (t(z ~ Az < cosd < z + Az) 
~z 2Az 



(14) 



2.3 Special final states 

The formation of x^H ions in reaction (|l^) is an arti- 
fact of our model due to the treatment of the electrons 
as fixed charge distributions. The cross-sections for these 
processes turn out to be small, and usually the final rif is 
small, so that the electron screening is not very important. 
For the purpose of cascade calculations, one can count the 
a;~H formation as the x~p events with the corresponding 



values of n/. If, cos 9, and ZiE'targct- Another channel in- 
volving a;~H ions is related to the formation of mctastable 
molecular states like 



p{x iVfnfl 



(15) 



where a deeply bound a;~H ion forms a loosely bound state 
with the proton. These molecular states can be rather 
stable and often do not dissociate within a reasonable 
amount of computer time. In our calculations we consider 
the metastable molecular states as final states. We used 
the following criteria for the metastability: first, the colli- 
sion time must exceed 



tmo\ = 50/ Winit 



(16) 



where Winit is the initial velocity of the x~p in the labora- 
tory system. With the choice of the time interval (|T6|), the 
colliding particles reach their asymptotically free final tra- 
jectories for most non-resonant collisions. Second, the x~ 
must form a bound state with one of the hydrogen atoms 
and the binding energy must not vary by more than 1% 
within the time 

r = 20^ (17) 

which corresponds to 20 classical periods of the initial x~p 
atom. Once metastability is reached, the event is counted 
as an a;~H event. 



3 Results 

The classical-trajectory Monte Carlo method described in 
Section ^ has been used to obtain the coUisional cross- 
sections needed in calculations of the cascades in fi~p, 
TT^p, K~p, and pp. The same method can also be used 
in a direct simulation of the atomic cascade without us- 
ing pre-calculated cross-sections. For fx~p and pp atoms 
experimental data at low density are available for direct 
comparison with the cascade calculations |l^. We will, 
therefore, present detailed results for these two cases. The 
initial stages also affect the cascades in tt~p and K~p be- 
cause they determine the kinetic energy distribution in 
the intermediate stage of the cascade where nuclear ab- 
sorption becomes important. 

The calculations have been done for Ui =8 — 20 for 
Hi = 13 — 35 for pp and 9 values of the laboratory 
kinetic energy in the interval 0.05 eV < T < 20 eV. At 
T = 1 eV the cross-sections have been calculated down 
to = 4 for and Ui = 8 for pp. For each initial 
state {ni,T), 1000 classical trajectories have been calcu- 
lated as described above. The orbital quantum number li 
was distributed according to the statistical weight. For the 
purpose of illustration, a larger number of trajectories (up 
to 10000) have been calculated for some initial states in 
order to reduce statistical errors. Preliminary results have 
been shown in ||ll| . 

We compare the results of the classical Monte Carlo 
(CMC) calculations with those of the semiclassical approx- 
imation. Bracci and Fiorentini [fl3| calculated the Coulomb 
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cross-sections for muonic hydrogen scattering from atomic 
hydrogen in a semiclassical model. Though the approach 
IP^ may be unsuitable for treating the low n states, where 
more elaborate calculations give much smaller values for 
the cross-sections [Q, it can be expected to give a fair 
description of the high n region. In the case of Stark mix- 
ing we use the fixed field model JT4| for comparison. In 
the case of molecular target, we obtained a semiclassical 
estimate of the Stark cross-sections by using the spher- 
ical symmetric electric field corresponding to the charge 
distribution of a II2 molecule in the ground state. 
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♦ CMC, H2 
□ CMC, H 

— Bracci, Fiorentini 




3.1 Muonic hydrogen 
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3.1.1 Coulomb deexcitation Fig. 5. The energy dependence of the Coulomb cross-sections 

for muonic hydrogen with rii = 13 and molecular (filled dia- 
The n dependences of the total cross-sections of the Coulomb monds) and atomic (squares)^ hydrogen target. The error bars 
deexcitation for collisions with molecular and atomic hy- 
drogen 



are statistical. The curve is the semiclassical result from |l^ . 



P)niU + H2 



fi-f 
If 



fi-f ' 



H + H 



with nj < Hi are shown in Figure ^. The cross-sections 
increase steadily with increasing n as the fj,~p becomes 
larger and the energy spacing between the n levels smaller. 
The cross-sections for the atomic target at the laboratory 
kinetic energy T — 1 eV are very close to the semiclassical 
results of Bracci and Fiorentini ||l^ . The cross-section for 
the molecular target is larger by a factor of about 2 — 3. 
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for the atomic targ et is in fair agreement with the semi- 
(18) classical result |12| for energies above 1 eV. The energy 
dependence of the CMC cross-sections is approximately 

(19) 

given by 1/VT corresponding to constant rates. This is in 
contrast to the 1/T behavior found for low energies in 

The distribution over final states nj is completely dif- 
ferent for the molecular and the atomic targets as illus- 
trated in Figure ^ showing the ^-average cross-sections 
for fipatl eV. The calculations for atomic target 
predict that An — 1 transitions dominate the Coulomb 
deexcitation in agreement with the semiclassical result 
For the molecular target, the transitions with An > 1 
are strongly enhanced as compared to the atomic case. 
The shape of the n/ distribution depends on the initial 
state nf. with decreasing it becomes narrower and its 
maximum shifts towards smaller values of An. For = 
13, the transitions An = 2 — 3 dominate. Figure |^ shows 
the nj dependence for initial state nt = 9: the transitions 
with An = 1 are most likely, but the An > 1 transitions 
still make up a substantial fraction of 38% of the Coulomb 
cross-section as compared to 19% for atomic target. 



3.1.2 Stark mixing and elastic scattering 

The Stark collisions change the orbital angular momentum 
while preserving the principal quantum number: 



Fig. 4. The n dependence of the muonic hydrogen Coulomb 
cross-sections at the laboratory energy T — 1 eV for molecular 
(filled diamonds) and atomic (squares) hydrogen target. The 
curve is the semiclassical result from |13]. 



An example of the energy dependence of the total 
Coulomb cross-sections {nf < ni) for = 13 is shown 
in Figure I The cross-sections calculated with molecular 
target are approximately twice as large as the atomic ones 
in the whole energy range considered. The CMC result 



p)yuU + H2 
{x^P)n,h + H 



{x p)r 
ix^P)n, 



kh; (H + H), 
H. 



(20) 
(21) 



The CMC results for the n dependence of the Z-average 
Stark mixing cross-section are shown in Figure |^. The 
Stark cross-sections calculated with molecular target are 
less than twice the atomic ones. This is due to two reasons. 
First, there is a considerable molecular screening effect 
because the electric fields from the two hydrogen atoms 
partly cancel each other. Second, the Coulomb cross-section 
makes up a larger fraction of the total cross-section in the 
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Fig. 6. The n/ dependence of the Coulomb cross-sections for 
muonic hydrogen with Ui = 13 and laboratory kinetic energy 
T = 1 eV for collisions with molecular (filled diamonds) and 
atomic (squares) hydrogen target. The semiclassical result from 
I p^if is shown with filled circles. 
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Fig. 7. The n / dependence of the Coulomb cross-sections for 
muonic hydrogen with n; = 9 and laboratory kinetic energy 
T = 1 eV for collisions with molecular (filled diamonds) and 
atomic (squares) hydrogen target. The semiclassical result from 
is shown with filled circles. 
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Fig. 8. Stark cross-sections for muonic hydrogen for molecu- 
lar (filled triangles) and atomic (circles) hydrogen target. The 
curves show the results of the fixed field model for molecular 
target (solid line) and atomic target (dashed line). The labo- 
ratory kinetic energy is T = 1 eV. 
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Fig. 9. The energy dependence of the Stark cross-sections for 
muonic hydrogen in the state n = 9 for molecular (filled tri- 
angles) and atomic (circles) hydrogen target. The curves show 
the results of the fixed field model for molecular target (solid 
line) and atomic target (dashed line). 



molecular case. The classical Monte Carlo results for the 
atomic target are in a good agreement with the semiclas- 
sical fixed field model. At low n, where the inelasticity due 
to the Coulomb deexcitation is small and can be neglected 
in the calculation of the Stark cross-sections, there is a 
good agreement between the classical Monte Carlo results 
for the molecular target and the corresponding semiclas- 
sical model. 

Figure ^ shows the energy dependence of the Stark 
cross-sections for n = 9. The classical-trajectory model 
and fixed field model are in agreement with each other 
for kinetic energies above 10 eV (molecular target) and 
2 eV (atomic target). At lower energies where the Coulomb 
transitions make up a substantial part of the cross-sections, 
the fixed field model overestimates the Stark cross-sections. 



The Stark mixing and elastic scattering processes, (|2C| ) 
and (|2l|), lead to a deceleration of the exotic atom. Their 
importance in the kinetics of atomic cascade can be esti- 
mated with the corresponding transport cross-section 

= |(l-cos0)^^dr2 (22) 

where d(T„_^„/di7 is the differential cross-section for the 
processes ( pO| ) or (^ averaged over I. This estimate based 
on the transport cross-section neglects the Coulomb deex- 
citation process which can lead to both deceleration and 
acceleration, and, in the case of molecular target, the ad- 
ditional deceleration due to excitation of the H2 molecule. 
The n dependence of the transport cross-sections at 1 eV 
for muonic hydrogen scattering from hydrogen atoms and 
molecules is shown in Figure y^. There is a fair agreement 
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Fig. 10. The n dependence of the transport cross-sections 
for muonic hydrogen at T = 1 eV. The results of the classical 
Monte Carlo model with molecular (filled circles) and atomic 
target (circles) are shown in comparison with the semiclassical 
fixed field model. 
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Fig. 11. Differential cross-section Aa/Az for muonic hydrogen 
for rii = 13 and laboratory kinetic energy T = 1 eV. The 
classical Monte Carlo results for atomic target are shown with 
filled diamonds and the curve corresponds to the semiclassical 
fixed field model for atomic target. 



between the CMC and the fixed field model for atomic 
target below n ^ 8. For higher n, the inelastic effects 
due to the Coulomb deexcitation process become impor- 
tant, and the fixed field model overestimates the trans- 
port cross-section. For molecular target, the discrepancy 
between the two models is larger because the Coulomb 
cross-section makes up a larger fraction of the total cross- 
section as compared to the CMC model with atomic tar- 
get (for n = 10 and T = 1 eV the fractions are ^ 0.24 for 
molecular target and ~ 0.11 for atomic target). 

Figure |ll| shows the /-averaged differential cross-section 
(using 20 equally spaced bins in z) summed over all the 
final channels for rii = 13 in the classical Monte Carlo 
model with atomic target. The cross-section is in good 
agreement with that of the semiclassical fixed field model. 
The pattern of maxima and minima in the semiclassi- 
cal differential cross-sections is a characteristic feature of 
quantum mechanical scattering, which, of course, cannot 
be reproduced in a classical model. 

The kinetic energy of the x~p in the final state is im- 
portant for detailed cascade calculations. Let T'j.-p, Th 
and be the CMS kinetic energies of the x~p, the H 
(for atomic target), and the H2 (for molecular target). The 
total kinetic energy is shared among the two {x~p and H) 
or three atoms {x~p and two hydrogen atoms): 
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Fig. 12. Distribution over the ii~p kinetic energy as a frac- 
tion of the total energy in the final state, T^-^/Ttot, for the 
Coulomb deexcitation of muonic hydrogen in the initial states 
Hi = 9, 13 at the laboratory kinetic energy T = 1 eV. The 
phase space distribution is shown for comparison. The vertical 
arrows indicate the final energies of the two-body final 
states fj,~p + H (at.) and jj,~p + H2 (mol.). 



kinematical boundaries: 



Jtnt — 



T - 



X p 



Th, atomic target 

+ ^-Etargct , molecular target 



(23) 



T - 
0< 1^ 



< 



2Mh 



2Mh 



(25) 



In the case of atomic target, the energy of the x pin CMS 
is fixed: 



T 



X p 



M^p + Mh 



(= 0.47 for fi-p) 



(24) 



where A/h and M^p are the masses of the hydrogen atom 
and the x^p atom, correspondingly. The case of molecular 
target corresponds to a three-body final state with the 



The upper boundary (0.64 for muonic hydrogen) is reached 
when the hydrogen molecule remains in its ground state 
corresponding effectively to a two-body {{x~p) + (H2)) fi- 
nal state. Figure |l^ shows the distributions in T^-p/Ttot 
for Coulomb deexcitations calculated in the classical Monte 
Carlo model for muonic hydrogen with Ui = 9, 13 and 
T = 1 eV. The approximation of effective two-body final 
states clearly fails, whereas the pure phase space distribu- 
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X p \ \ 

gives a fair description of the results 
3.2 Antiprotonic hydrogen 



The atomic cascade in antiprotonic hydrogen starts around 
Hi ~ 30; thus classical mechanics is even a better ap- 
proximation than in the muonic hydrogen case. The n 
dependence of the Stark mixing, Coulomb deexcitation, 
transport, and the pR formation cross-sections is shown 
in Figure [l3|, and the energy dependence is demonstrated 
in Figure p^. As with muonic hydrogen, the fixed field 
model overestimates the Stark mixing and especially the 
transport cross-section because the inelasticity effects due 
to Coulomb deexcitation are not included in this frame- 
work. 
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Fig. 13. The n dependence of the cross-sections of the cascade 
processes in antiprotonic hydrogen calculated in the classical- 
trajectory Monte Carlo model for molecular target. The lab- 
oratory energy is T = 1 eV. The Stark mixing and transport 
cross-sections calculated in the fixed field model for molecular 
target are shown with solid and dashed lines, respectively. 



Figure shows the distribution over the final states 
n/ for the Coulomb deexcitation of the antiprotonic hy- 
drogen at the laboratory energy T = 1 eV. For high n 
initial states, the most probable Coulomb transitions are 
the ones with a large change of the principal quantum 
number (Z\n ^ 1), with the molecular target being es- 
sential for this feature. A very important consequence of 
this result is that at the beginning of the atomic cascade a 
small number of Coulomb transitions is sufficient to bring 
the pp to the middle stage, where, depending on the target 
density, the radiative or Auger deexcitation takes over. 

The dependence of the Coulomb cross-sections on the 
angular momentum li of the initial state is weak, see Fig- 
ure ^ for antiprotonic hydrogen with rii ~ 25. The Stark 
cross-sections show a moderate dependence on If. they are 
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Fig. 14. The energy dependence of the cross-sections of the 
cascade processes in antiprotonic hydrogen in the state rii — 25 
calculated in the classical-trajectory Monte Carlo model for 
molecular target. 
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Fig. 15. The n/ dependence of the Coulomb cross-sections 
for antiprotonic hydrogen for = 20, 25, 30 and the laboratory 
kinetic energy T = 1 eV. For the sake of clarity we do not show 
the statistical error bars. 



smaller for the circular states {li = — 1) than for the low 
li, by about 50%. The reason for this is that the elongated 
ellipses in the low I states are more easily perturbed by 
the electric field of the target molecule. A similar effect is 
expected if a quantum mechanical description of the pp is 
used: the size of the pp as estimated by the expectation 



value of r is given by 



(r^ 



(5n2 + l-3?(/ + l)) . 



(27) 



For high n states, the expectation value of for the cir- 
cular state is only 40% of that of the ns states. 
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Fig. 16. Tlie h dependence of the Coulomb (filled diamonds) 
and Stark (filled triangles) cross-sections for antiprotonic hy- 
drogen for rii = 25 and T = 1 eV. The results are calculated 
in the classical Monte Carlo model with molecular target. 



4 External Auger effect in the eikonal 
approximation 

In our treatment of the Coulomb and Stark mixing colli- 
sions in Section ^ the electronic degrees of freedom were 
assumed to be frozen. These degrees of freedom, however, 
play an important role in the Auger deexcitation process 



(28) 



The Auger transitions are often treated in the Born ap- 
proximation [0 that gives (conveniently) energy indepen- 
dent rates. However, this approximation violates unitarity 
for some important ranges of principal quantum numbers 
and kinetic energies. For kinetic energies in the range of 
few eV, the eikonal approximation |15|| provides a more 
suitable framework. In this section, we use the eikonal 
approach to calculate Stark mixing and Auger deexcita- 
tion simultaneously. As a result, the corresponding partial 
wave cross-sections are consistent with unitarity. 

The cross-section for the process ( p8| ) was calculated 
in |l5[| by assuming that the exotic atom moves along a 
straight line trajectory with constant velocity v through 
the electric field of the hydrogen atom at rest. The cross- 
section is given by 



Auger n 



Pip)pdp 



(29) 



where P{p) is the reaction probability for the impact pa- 
rameter p: 



Pip) = 1 



-Hp) 



V J-oo 



(30) 

The reaction rate, Pmii (R), at distance R is the sum of the 
partial rates Fnji^nfifiR) over all final states. According 
to Uai the estimated rates are 



Rmli^njlj [R) — 7 



1 



(il'2 + fe2)3 



-71- 



■exp(-2i?) (31) 



where fcg is the electron momentum, h = 1.5, and the 
parameters 7 and 71 are given by 



exp {{—A/kg) arctan k^,) 
3 (l + A:2)6(i„exp(-2Vfce)) 



7 = -^f^^ 



(32) 
(33) 



where C^^qiq is a Clebsch-Gordan coefficient and R^^i^ 
is the radial matrix element [p^ . The transition rate is 
proportional to the square of the dipole matrix element, 
therefore only transitions with Al — \l f — li\ = 1 are pos- 
sible. 

The Auger deexcitation rate, as a function of n, peaks 
at the n- value where the energy released in a An = 1 tran- 
sition is just sufficient to ionize the hydrogen atom. The 
effect of these high-rate Auger transitions is that the in- 
elastic cross-sections for some partial waves are not small 
in comparison with the unitarity limit. Therefore the cor- 
responding inelasticity should be taken into account in the 
calculations of other coUisional processes. One can expect 
that taking the Auger effect into account will reduce the 
other inelastic cross-sections. In order to examine this ef- 
fect, we include the Auger deexcitation in the framework 
presented in ||lj] for calculating Stark mixing and elas- 
tic scattering. In the same way as the nuclear absorption 
processes in hadronic atoms were taken into account via 
imaginary energy shifts of the ns-states, the Auger deex- 
citation process is included via the imaginary absorption 
potential, — zr'„;(i?)/2. The calculations can be done in 
the close-coupling model, the semiclassical model, and the 
fixed field model. In the case of the fixed field model, the 
time-dependent Schrodinger equation for the set of the 
linear independent solutions forming the x matrix 
A is given by 

iA{t) = H{t)A{t) (34) 
where the interaction is given by 



H{t) = Z 



1 



R^{t) 



AE 



(l + 2i?(i) + 2ii'2(t))e-2«(*) 

r{R) 



(35) 



Here AE is a diagonal matrix corresponding to the en- 
ergy shifts due to the vacuum polarization and the strong 
interaction. The term r{R) is a diagonal matrix with the 
matrix elements Fni {R) . The factor Z originates from the 
dipole operator and has the following matrix elements 
{t=\nlA),j^\n{l-l)A)): 



Zij — 



' {P - A^){n^ - P) 



3n 



2p^p\ (2/ + l)(2Z-l) 



(36) 



The solution of Equation ( |34| ) using the method de- 
scribed in gives the scattering matrix 5^^. The cross- 
sections for the transitions nJi riilf are given by 



1 



2k 



1 fc2 ^ 
,7 



(2J+ 1) 
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Fig. 17. The energy dependence of the Auger deexcitation 
rates for muonic hydrogen in liquid hydrogen. The results of 
the eikonal approximation are shown with solid lines and those 
of the Born approximation with dashed lines. 




Fig. 18. The n dependence of the Auger deexcitation and 
Stark mixing rates at 1 eV for muonic hydrogen in liquid hydro- 
gen. The results of the eikonal multichannel model are shown 
with solid lines. The Auger deexcitation rates calculated in the 
the Born approximation and the Stark mixing rates obtained 
without Auger deexcitation are shown with dashed lines. 
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J 

Fig. 19. The J dependence of the /-average partial wave cross- 
sections for muonic hydrogen for n = 7 and laboratory kinetic 
energy T — 3 eV. The cross-sections for Auger deexcitation 
with An — 1 are shown with diamonds, those of Stark mixing 
with histograms, and the unitarity limit with a thick solid line. 
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Fig. 20. The initial / dependence of the Auger deexcitation 
cross-sections at 1 eV for muonic hydrogen at n = 7. The result 
of the eikonal multichannel model is shown with filled circles 
and that of method with filled squares. 



X \{nf,JMAlf\S'''' ~l\nf,JMAk)\''{37) 
and the ones of the Auger deexcitation by 



Auger 



21 



(2j + i)((2;, + i) 

n, ■ J MAI flS^^l n, ; J MAI ^ ' ^ 



1 ^ 

J 



Alf 



We will refer to this framework as the eikonal muhichannel 
model. 

The Z-average Auger deexcitation cross-sections calcu- 



lated with the method of 
agree closely with our resu 



_ (Equations (|29|) and (^) 
ts in the eikonal multichannel 
model. Figure |l^ shows the /-average An = 1 Auger deex- 
citation rates in muonic hydrogen for n = 3, 5, 7. The rates 
have been calculated in the eikonal approximation and the 



Born approximation. The rates in the eikonal approxima- 
tion are lower in the low energy range, but they approach 
the ones of the Born approximation for high energies. The 
n dependence of the Auger deexcitation and Stark mixing 
rates for muonic hydrogen is presented in Figure The 
two approaches are in a fair agreement with each other ex- 
. cept for the states n — 6,7 where the Auger rates have the 

) I j . (38) highest values. For the state n = 7, the Stark mixing rates 
are reduced by almost 50% when the inelasticity due to the 
Auger effect is included. This resembles the situation with 
the eikonal and the Born approximations which disagree 
when the Auger deexcitation cross-sections are large, in 
which case the eikonal approximation gives smaller cross- 
sections than the Born approximation. The explanation 
of this effect is given in Figure ^ showing the average 
partial wave cross-sections for the collision {fip)n=7 + H. 
The Auger deexcitation cross-sections are saturated in the 
low angular momentum region and, therefore, the Born 
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approximation fails. Though the Z-average results agree 
for the two eikonal approaches, the I dependence of the 
cross-sections in the eikonal multichannel model is weaker 
because of the effect of Stark mixing as demonstrated in 
Figure |0[ 

The eikonal approximation as described above does not 
give the differential cross-section and distribution over fi- 
nal // for the Auger transitions. The partial wave cross- 
sections. Figure show that the main contribution to 
the Auger cross-section comes from low partial waves, i.e. 
from the strong mixing region. This suggests that the dis- 
tribution iu-lf is nearly statistical and that the differential 
cross-section is less forward-peaked than the elastic and 
Stark differential cross-sections Q . 

5 Conclusions 

The collisional deexcitation mechanisms of the exotic hy- 
drogen atoms in highly excited states have been investi- 
gated in detail using the classical-trajectory Monte Carlo 
method. The Coulomb transitions have been shown to 
be the dominant mechanism of collisional deexcitation of 
highly excited exotic atoms. Target molecular structure 
has large effects on the Coulomb deexcitation. In par- 
ticular, the distribution over the final states favors large 
change of the principal quantum number n contrary to the 
case of atomic target. This feature is very important for 
the cascade kinetics as it leads to a fast deexcitation and 
a significant acceleration at the initial stage of the atomic 
cascade . The calculated cross-sections provide a more 
reliable theoretical input for further cascade studies by re- 
moving the long standing puzzle of the so-called chemical 
deexcitation ||l|], which was used, on purely phenomenolog- 
ical grounds, in many cascade calculations without clari- 
fication of underlying dynamics. 

The external Auger effect has been studied in an eikonal 
multichannel model which allows us to calculate Stark 
mixing, elastic scattering, and Auger deexcitation simul- 
taneously. Partial wave cross-sections computed in this 
framework are consistent with unitarity. For ranges of 
principal quantum numbers and kinetic energies where the 
unitarity constraint is important, the Auger cross-sections 
computed in this model arc significantly lower than those 
of the Born approximation [Q| . 

The first results of cascade calculations using the cross- 
sections of and the present paper have been presented 
in [0|l^ . More detailed results of the cascade calculations 
will be discussed in a separate publication . 
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